s  iV>' 

■:iS  'V”  *£  » 


assrafp 

|i',,  f  1  ‘  I  ■> 


feup . 

■  6  ,4 

Stiff  • 

EwAv.k 

mkp  ,  &•*> 


i»w? 

$**'*•“  #V «S  1  Vv, 

-  '>£■* 


‘‘'■if.y 


'% 


rT&l  Kr  ‘  •  * 

"  V, 
:V  ,a£li 


mmi 


r  ? 


PI  |  m 
Pm  m$h  4  i  $ 

Spm 

PiMp« 

Mm®*  m  ■■  ■ 


*  '£yi'i. 
'■i  *  v 


uNP®  t  3Av  ffi 


;  '  1  _ 


r> 


APPROVED  POR  PUM  ^ASt* 
OISTMMiriON  UNLIMITED 


»d  p> 


UNCIASSIFIED 


DOCUMENT  CONTROL  DATA  R  &  D 

.*»«•:  liri/i'  (faisiWn  .ifiori  of  title,  /><>(/>  at  iihttrat  I  mid  index int*  .innotatmn  mu -t  he  enteted  wfirn  f/i«*  i/vrruf/  report  is  c  las  si  tied) 


\  ONlbtNA  IINC  activity  f  Ciifp0M/p  »n  ithor)  30,  REPORT  Sf  C  JRITy  C  L  A  SSI  f  I  r  A  T  ION 

UNCIASSIFIED 

Naval  Ordnance  Laboratory  |77>  GROUP 


Naval  Ordnance  Laboratory 

White  Oak,  Silver  Spring,  Maryland  20910 


HM-ORT  T  I  T  L  fr 


Finite  Amplitude  Propagation  of  an  Underwater  Explosion  Shock  Wave  Along  a 
Strongly  Refracted  Ray  Tube 


4  OtSCRiP  nvf  N  O  T  F  s  (  Type  of  rep  >rt  and  nidus  i  ve  dates  j 


au  TmoRiSI  ff'ir.f  name*,  middle  mitiml,  taxi  name) 

John  F.  Goertner 


fe  Rfc  PO«  T  DATE 

1  December  1971 

ra.  TOTAL  NO  op  pacts 

ill  ♦  24 

?b  NO  OT  RFFS 

4 

0«  CONTRACT  ON  GRANT  NO 

9 a.  ORlG!NATOR*5  REPORT  NuMR£RlS|  ] 

h  PHOjrC  ’  NO 

DRA  Subtask  NA  002/20 

c 

NOL  TR  71-139 

•it  other  report  no<S>  (Any  other  number*  that  may  be  assigned 

d 

fhi*  report > 

Dl  5  1  Rt  0J  T  IQI,  STATEMEN1 


Approved  for  Public  Release;  Distribution  Unlimited 


*3  SPOSSCNiNG  Mlli’AWi  ACTIVITY 

Defense  Nuclear  Agency 

Washington,  D.C. 


The  effect  of  wave  overtaking  on  the  peak  pressure  of  refracted  underwater  explosion 
shock  waves  is  estiaated  by  means  of  a  finite  amplitude  calculation  for  the  wave 
propagation  along  refracted  ray  tubes  given  by  simple  ray  theory.  The  results 
Indicate  that  a  reduction  in  peak  pressure  occurs  due  to  increased  wave-overtaking 
along  refracted  ray  pathe  to  first  caustics.  Reductions  In  the  peak  pressure  no 
greater  than  5 i  occurred  for  the  particular  refractive  conditions  treated,  but 
preliminary  work  with  other  conditions  indicates  much  greater  reductions  can  occur. 


DD  /r.,1473 


S'S  010  1. 807. 6  80  t 


(PAGE  1) 


UWCIASSiriED 

Security  Clessif  teal  ton 


HOLE  WT 


Underwater  Explosions 

Refraction 

Shock  Wave 

Finite  Anplltude  Wave  Propagation 
Caustics 

Nuclear  Explosions 
Ray  Tracing 


DD  :r..1473  back, 


(PA(,r  m 


UNCLASSIFIED 

Security  Claiuflcation 


NOLTR  71-139 


FINITE  AMPLITUDE  PROPAGATION 
OP  AN 

UNDERWATER  EXPLOSION  SRCCK  WAVE 
ALONG  A  STRONGLY  REFRACTED  RAY  TUBE 


Prepared  by: 
John  F.  Coertner 


ABSTRACT:  The  effect  of  ware  overtaking  on  the  peak  pressure  of  refracted 
underwater  explosion  shock  wares  Is  estlaated  by  aeans  of  a  finite  amplitude 
calculation  for  the  ware  propagation  along  refracted  ray  tubes  given  by  slnple 
ray  theory.  The  results  indicate  that  a  reduction  In  peak  pressure  occurs  due  to 
Increased  wave-overtaking  along  refracted  ray  paths  to  first  caustics. 

Reductions  In  the  peak  pressure  no  greater  than  5%  occurred  Tor  the  particular 
refractive  conditions  treated ,  but  prollalnary  work  with  other  conditions 
indicates  such  greater  reductions  can  occur. 
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Finite  Amplitude  Propagation  of  an  Underwater  Explosion  Shock  Wave  Along  a  Strongly 
Refracted  Ray  Tube 


Present  understanding  of  the  refraction  of  underwater  explosion  shock  waves  has 
been  achieved  largely  through  the  use  of  acoustic  theory,  i.e.,  theory  of  low 
amplitude  pressure  waves,  with  an  overall  modification  to  account  for  explosion 
pulse  propagation.  The  work  reported  here  Is  the  first  part  of  a  study 
undertaken  to  determine  whether  it  is  necessary  to  modify  the  present  methods  so 
as  to  intrinsically  Incorporate  the  non-linear  effects  due  to  the  finite  amplitude 
of  explosion  pulses.  The  work  was  supported  by  Defense  Atomic  Support  Agency 
(now  Defense  Nuclear  Agency)  Subtask  RA  002/20,  Underwater  Shock  Theory  /Energy 
Focusing  and  Refraction  Effects. 
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the  course  cf  this  work. 
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1.  DWRODOCFIOW 

A  general  nathenatloal  solution — either  analytic  or  numerical— to  the  pr atlas 
of  refraction  of  an  underwater  explosion  shock  were  by  a  strong  velocity  gradient 
does  not  yet  exlet.  Today's  solutions  consist  of  Judicious  patching  together  of 
various  asymptotic  solutions  and  experlaental  results.  For  exmspie,  Brockhurst, 
huN,  end  Arons  (1961)*  demonstrated  the  usefulness  of  acoustic  ray  thtcary  In 
describing  refracted  vers  pressures  there  strong  focusing  has  not  occurred.  Near 
regions  of  strong  focusing  In  the  ocean  (convergence  sans  caustics),  Blatsteln 
(1971)  has  demonstrated  the  usefulness  of  a  modified  wave-length-dependent  ray 
thee sy  to  describe  the  experimentally  measured  pressure  signatures.  However, 
neither  simple  ray  theory  nor  the  modified  ray  theory  take  Into  aceoost  the  wavs 
overtaking  effects  Inherent  to  shock  wave  propagation.**  Tbs  above  mentioned 
authors  in  applying  their  respective  treatments  Incorporated  the  Judicious 
assmqrtlon  that  the  effects  of  wavs  overtaking  are  the  same  as  far  a  non-refracted 
underwater  explosion  shock  wavs  that  has  traveled  the  same  distance.  Tbs  purpose 
of  this  note  1*  to  check  this  assumption  by  calculating  the  propagation  of  an 
underwater  explosion  shock  wave  along  typical  acoustically  refracted  ray  tubes. 

To  do  this  we  use  a  method -of -characteristics  niasrloal  Integration  of  the  fluid 
dynamical  equations  expressing  conservation  of  ness  and  of  momenta*  In  the 
direction  of  the  ray  tube. 

This  paper  consists  of  three  sections.  The  first  describes  the  approximate 
thefgy  end  lists  the  pertinent  equations.  Tbe  second  gives  a  brief  account  of 
the  essential  features  of  the  moerlcal  calculations.  And,  the  third  soonrlsss 
the  results  calculated  for  a  sampling  of  refracted  ray  tubes  for  a  source  depth 
and  velocity  profile  yielding  a  thermooline-related  caustic. 

2.  THEORY 


The  problem  considered  Is  the  lmrlscld  unsteady  compressible  flow  along  the 
axis  of  a  rigid  tube  of  arbitrarily  changing  area. 

Constitutive  Equations*  Tbe  equations  used  to  describe  the  notion  are  the 
roieirtun  ana  contlmiuty  equations  and  a  pressure-density  relation  for  the  water. 


•Refer  to  list  of  references  on  Page  2k 

**Fo*  discussion  of  wave-overtaking  in  underwater  explosion  shock  wave  propagation 
see  Snay,  1966,  Section  IF  or  Cole,  19WJ,  Sections  2. 1-2. 6 
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‘UUag  the  axis  of  the  ray  tuba  tha  ■osnntvai  aquation  radueaa  to 

|H  ♦  U|S  ♦  i|B  .  0 

at  ox  ^  bx 

and  tha  continuity  aquation  baconaa 


t  ■  tine 

x  ■  diatanea  along  tha  axia 
u  ■  part  ic la  velocity 

p  -  praaaura  (relative  to  tha  aahlant  value) 
f*  ■  uatar  density 

A  ■  A(x),  araa  of  tha  ray  tuba  cross-section 

Par  tha  pressure-density  ralatlon  wa  uaad  tha  nodifiad  Tait  aquation 
(Cola,  19W,  page  U»). 

P  - 


Utara  e  ■  ambient  sound  speed 

0  (subaerlpt  "o"  lndloataa  anblant  state) 
n  •  constant 

With  n»6,0  aquation  (3)  dascrlbas  tha  aechanlcal  behavior  of  tha  uatar  behind  both 
shocks  and  laantropas  to  within  about  *  0.50  for  prassuraa  in  tha  range,  0sp<10,000 
pal;  i.a. ,  flow  par  east  era  such  as  tha  density,  sound  speed,  particle  velocity, 
Aock  velocity,  and  Rleaann  function  are  a st lasted  to  be  within  about  to.  50  of 
their  true  values.* 

Essentially,  by  Introducing  aquation  (3)  to  describe  tha  entire  flow,  us  a re 
neglecting  tha  effects  of  tha  snail  but  variable  entropy  lncreaaes  iftlch  occur  as 
tha  water  Is  traversed  by  tha  outwardly  propagating  riiock  front.  In  this  report 
we  will  consistently  asks  use  of  tha  simplifications  ifelch  follow  fron  this 
assiaptlon  of  unlfara  entropy  throughout  tha  entire  flow.  Aa  pointed  out  above, 
this  leads  to  a  at  lasted  errors  of  tha  order  of  10.50  In  tha  computed  flow 
paraastara. 


icterlstlc  Equations. 


Introducing  tha  Rleaann  function 
1 


•  7  h  * 


•These  accuracy  Units  were  estimated  fron  Tables  3*3  and  3.*  of  NOME  70-31  (H.  0.  Snay 
and  A.  R.  Krlebel,  1970)  £y  making  use  of  "Zeta  equation"  as  an  intermediate  standard. 
Using  n-6.0,  equation  (3)  is  essentially  equivalent  to  the  Zeta  equation;  and,  over  the 
pressure  range,  zero  to  10,000  pel,  the  calculated  values  behind  shocks  and 
lsentropes  far  all  the  aechanlcal  flow  parameters  are  the  sane  to  within  *0.20. 
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and  noting  that 


equations  1  and  2  can  be  transformed  to 

^  *  («*) 

^  *  («) 

In  deriving  equations  (6)  and  (7)  ve  have  implicity  used  the  uniform  entropy 
assumption. * 

Equations  (6)  and  (7)  are  the  so-called  “characteristic  equations."  They  give 
the  time  rates  of  change  of  the  quantities  (u+o)  and  (u-o)  to  observers  traveling 
in  the  x-t  plane  at  velocities  (u+c )  and  (u-c),  respectively.  Given  suitable 
boundary  conditions  they  can  be  integrated  numerically  throughout  the  region  in  the 
x-t  plane  between  the  outward- propagating  shock  front  and  some  starting  point 
(value  of  x)  along  the  ray  tube. 

To  solve  equations  (6)  and  (7)  and  also  to  calculate  the  pressure  we  need 
several  additional  relations.  The  following  are  derived  from  equations  (3)/  00, 
and  (5): 


=  & 


ox 


uc  BA 
A  Bx 


d(u-o) 

Bx 


uc  BA 
A  Bx 


(5) 


(6) 

(7) 


Hound  speed 

Pressure 

Riemsna  function 

Shock  front  boundary  condition 

Shock  front  velocity 


(8) 

(9) 

(10) 
(U) 
(12) 


^far  elucidation  of  the  role  of  entropy  and  the  significance  of  the  Riemann 
function,  the  reader  is  referred  to  NOLTR  65-52,  H.  G.  Snay,  1966,  pages  h-6-50. 
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Equation  (ll)  follows  from  (7)  neglecting  the  entropy  change  across  the  shock. 
Equation  (12)  Is  derived  from  the  mechanical  shock  relation,  p=()QuS,  (conservation 

of  ornss  and  momentum)  using  equations  (9)  and  (ll)  and  expanding  the  term 
2n 

n-l  /c-c  \ 

in  (9)  as  a  function  of  ^  ~c'  "J  811(1  dropping  terms  of  3rd  and  higher 

order. 

3-  lfUMERICAL  CALCUIATIOHS 


For  numerical  integration  equations  (6)  and  (7)  were  approximated  by  the 
following  difference  equations: 


A+  (u+ct)  ■  - 


tT  (u-o)  ■  + 


uc 


uc 


dlnA 

dx 


dlnA 

dx 


At 


ave 


At 


ave 


(13) 

(lfc) 


where 

A  *  change  between  mesh  points  along  a  positive  (u+a) 
character  1  st  1c . 

A"  »  change  between  mesh  points  along  a  negative  (u-a) 
characteristic. 


"Average"  refers  to  the  end  (or  mesh)  point  values.  The  characteristic  mesh  is 
approximated  by  line  segments  of  slope  ( u+c )aye  and  (u-c)ay#. 

For  a  more  complete  description  of  the  method  of  computation,  the  reader  Is 
referred  to  Goertner,  196%  Appendix  B.  In  this  section  we  will  discuss  details 
\dilch  are  unique  to  these  computations. 


Initial  and  Boundary  Conditions.  To  start  the  computation  (initial  condition) 
and  for  the  boundary  condition  at  the  starting  ray  tube  section  we  used  the  shock 
wave  similitude  equations  derived  from  experimental  data  for  TlfF  and  Pentollte 

(Cole,  1948,  page 8  238-242): 


rmax 


(15) 


0.06O 


(wl/3/s) 


-0.18 


(16) 


where  *  Initial  shock  pressure  In  pounds  per  square  Inch 

W  =  weight  of  TNT  in  pounds 
R  =  range  In  feet 

6  -  time  constant  of  the  exponential  pressure  pulse  In 

milliseconds. 
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k  -  2.16  x  1 (f  tar  TUT 
■  2.25  x  1 <f  tar  Fsntollt# 

Figure  1  illustrates  how  this  vu  dona.  Using  aquation  (15 ),  the  range 
corresponding  to  PMXa10,000  pal  la  calculated.  In  the  figure,  xQ  Is  the  starting 

point  for  the  computation;  this  la  set  equal  to  the  range  corresponding  to  p  ■ 

10,000  pel.  At  x*xq  the  pressure,  and  from  this,  the  Rleuann  function.  Is 

calculated  uelng  the  tine  constant,  ft,  given  by  equation  (l 6).  This  Is  the 

boundary  condition  used  at  x»x  . 

o 


To  start  the  oonputatlon,  I  suitably  spaced  nesh  points  are  placed  ala  (  the 
shock  path  through  the  point  (xQ,  0)  corresponding  to  the  slallarlty  equation  (15). 

For  these  computations,  the  If  Initial  mesh  points  were  Inserted  along  the  shock  in 
the  Interval  from  x  ■  xQ  to  1.1  xQ.  For  simplicity,  the  sketch  Illustrates  only 

4  Initial  mesh  points;  generally  a  greater  maker  -  32,  64,  or  128  -  were  required. 
Using  equations  (9)  through  (12)  we  can  calculate  the  shock  path  and  all  pertinent 
flow  variables  along  the  starting  shock. "  Starting  at  the  first  inserted  mesh 
point,  the  masrloal  integration  then  proceeds  from  shock  front  to  x*x  along 
successive  negative  characteristics. 


Comveraence  of  the  numerical  Computations.  In  the  course  of  development 
check-out,  three  versions  of  the  computer  program  were  used.  The  essential 


and 


difference  among  these  was  In  how  the  tern 


in  equations  (13)  and  (14) 


ave 


computed  ( approximated )  and  In  the  degree  of  control  over  this  approximation. 


Version  1  was  specifically  written  to  calculate  the  shock  wave  propagation 
along  a  spherically  diverging  ray  (l.e. ,  a  cone).  For  this  case  din  A/dx-2/x  so 
that  In  solving  equations  (13)  and  (14)  the  term  in  brackets  could  be  evaluated  at 
the  appropriate  mesh  points  and  then  averaged. 

Version  2  was  written  for  the  ray  tube  of  arbitrarily  changing  area;  A(x) 
being  calculated  by  3-polnt  Ugranglan  Interpolation  of  tabulated  values.  This  led 
to  a  slightly  different  approximation,  namely. 


uc 


(17) 


vdiers  the  "bar"  Indicates  average  and  "A"  the  difference  of  the  appropriate  mesh 
point  values. 


The  result  of  these  different  representations  of 


ainA 

dx 


In  equations  (13) 
ave 


and  (14)  on  the  shock  front  pressure  calculated  far  a  spherically  spreading  shock 
Is  shown  In  Figure  2.  We  see  that  as  If  increases  the  pressure-vs-dlstance  curves 
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INITIAL  PRESSURE  (PSI) 
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DISTANCE,  x  (FT) 

FIG. 2  VARIATION  OF  CALCULATED  PRESSURE-DISTANCE  WITH  MESH  SIZE 
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calculated  using  these  different  approximations  conrerge  fro*  opposite  sides  to 
what  Is  apparently  the  true  solution  to  equations  (6)  and  (7).  Thus,  by  good 
fortune,  us  obtained  a  valuable  check  on  the  computations. 

Control  of  the  Divergence  Term.  As  the  computation  proceeds  the  shock  front 
moves  outward  and  the  separation  of  adjacent  mesh  points  Increases  (Figure  l). 

As  a  result,  the  approximation  to  ^  A  computed  from  parameters  at  the  adjacent 


mesh  points  becomes  poorer  and  the  numerical  solution  degenerates  as  can  be  seen 
In  Figure  2.  Oxer  Initial  solution  to  this  problem  was  simply  to  Increase— double, 
re-double,  etc.— N,  the  number  of  Initial  negative  characteristics,  until  the 
numerical  solution  had  apparently  converged.  This  worked— except  for  some  ray  tubes 
which  diverged  rapidly  upon  entering  a  partial  shadow  zone.  For  these  the  rate  of 
convergence  with  Increasing  N  was  too  slow  and  computer  cost  became  excessive. 


To  compute  along  these  rays  we  designed  Version  3  which  subdivided  the 
characteristic  net— Introduced  additional  characteristics  and  discarded  others— as 
the  computation  proceeded.  The  modifications  which  constitute  Version  3  stem  from 
the  following  observations: 

(1)  That  the  computation  time  Is  roughly  proportional  to  the  mrniber  of  mesh  points 
calculated. 

(2)  That  In  both  Versions  1  and  2  the  major  source  of  systematic  computational 
error  appears  to  be  the  finite  difference  approximation  to  din  x  (Equations 
13  and  lk). 


Thus,  In  Version  3  we  attempted  to  calculate  the  minimum  number  of  mesh  points 
consistent  with  some  preaoslgned  level  of  overall  computation  error.  We  did  this 
by  monitoring  Ain  A  from  point  to  point  along  the  shock  trajectory.  In  the  event 

Alnftl  >  JaIuaJ  (IB) 


xtiere  pMnAJ^^  is  some  preset  value,  the  computation  switched  to  a  subroutine 

xdxich  introduced  N/2  Intermediate  mesh  points  along  the  half  of  the  previous 
negative  characteristic  adjacent  to  the  shock,  discarded  the  W/2  mesh  points  from 
the  other  half,  did  the  necessary  bookkeeping,  and  then  returned  to  re-calculate 
the  next  point  on  the  shock  trajectory.  If  necessary,  this  procedure  was  repeated 


until 


AlnA 


(19) 


Figure  2  Includes  a  segment  of  a  pressure-vs-dlstance  curve  calculated  using 
Version  3  with  LMnAJ  „  set  equal  to  0.2;  It  falls  between  the  curves  calculated 
using  Versions  1  and  ;  TmxWith  version  3  the  shock  trajectory  is  essentially 
independent  of  N,  the  number  of  Initial  characteristics.  (Using  this  program,  the 
value  specified  for  H  essentially  determines  the  area  of  integration  In  the  x-t 
plane,  and  Is  generally  specified  In  accordance  with  the  required  p-t  signature 
duration  at  some  value  of  x. ) 


The  value, 


AlnA 


>0.2,  was  selected  for  the  ray  tube  computations 

max 


summarized  in  Part  4.  This  value  was  selected  by  Inspection  of  the  csilculated 
pressure-vs-dlstance  data  such  as  shown  In  Figure  3  and  Table  1. 
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DISTANCE,  x  (FEET) 


FIG.  3  VARIATION  OF  CALCULATED  PRESSURE-DISTANCE  WITH  REZONING 
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TAME  1 


CHAHQE  IK  CALCUIATED  SHOCK  FROlfT  PRESSURE,  p  , 


as  a  rowcrioii  or 


x  *  200  ft  (vmx  ■  216  psl ) 

8.38  lbs  pentollts 
32  Initial  characteristics 
spherical  syatry 


AlnA  I  Relative  Error 


_  rsax 

(i) 

0.1*  — 

0.2  *0.3 

0.4  +1.5 

0.8  +4.6 

(Hot  Controlled )  +13.6 


•for  this  case  6k  initial  characteristics  sere  required  in  order  to  confute  out 
to  x  -  200  feet. 
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4.  RlSUIffib  FOR  REFRACTED  RAY  TUBES 


We  started  this  study  by  looking  at  refracted  rays  corresponding  to  the 
flooded  quarry  experiment  of  Brockhurst,  Bruce,  and  Arons  (l96l).  Figure  4  shows 
the  measured  average  sound  velocity  profile  and  the  corresponding  ray  diagram 
computed  by  ray  theory.  The  locations  of  the  rays  used  are  indicated  by  the  symbols, 
"A"  through  "E.” 

For  each  of  these  chosen  rays  ve  also  calculated  the  trajectories  of  two 
adjacent  rays  \hlch  defined  the  ray  tube  for  our  calculations.  Figures  5  through  9 
show  the  ray  tube  areas  for  these  rays  as  a  function  of  distance  along  the  central 
ray;  and,  compares  this  with  the  corresponding  area  for  a  non-refracted  ray  tube. 

Note  that  in  Figures  5  through  9  it  is  the  logarithms  of  the  ray  tube  areas  and 
path  lengths  tdilch  are  plotted.  This  is  appropriate  since  it  is  the  rate  of  change 
of  InA  (equations  13  and  14J  \&lch  governs  the  wave  propagation  along  the  tube; 
and,  it  is  the  ratio,  R/wA  ,  which  enters  into  the  boundary  and  initial  conditions 
of  the  calculation  (equations  15  and  l6).  This  is  also  in  accordance  with  the 
principle  of  similarity  for  underwater  explosion  shock  waves  (Cole,  1948,  page  110 ) 
tfclch  holds  for  spherical  waves  but  not  for  waves  propagating  along  refracted  ray 
tubes  (since  each  refracted  ray  tube  has  a  characteristic  length  of  its  own  which 
is  Independent  of  the  charge  size). 

We  now  make  use  of  the  shock  wave  similarity  principle  far  spherical  waves  to 
compare  the  ray  tube  area  functions  of  Figures  5  through  9.  We  do  this  by 
superposing  the  plots.  Figures  5  through  9,  and  sliding  them  along  the  straight 
line  representing  the  spherical  wave  until  ve  have  aligned  the  points  of  departure 
from  spherical  spreading.  (Note,  if  two  area  functions  can  be  made  to  coincide  in 
this  manner,  then  they  are  similar  in  the  sense  .that  a  calculation  far  one  can  be 
scaled  to  the  other— corresponding  values  of  wV3#  path  lengths,  and  times  will  be 
In  the  ratio  of  the  path  lengths  for  corresponding  points  of  the  area  functions). 

The  result  of  such  a  superposition  of  the  area  functions  is  shown  in  Figure  10, 
which,  far  the  purpose  of  this  report,  shows  the  essential  differences  among  these 
functions. 

Xir  selection  of  area  functions  (Figure  10)  is  characterized  by  two  extremes, 
case  "b"  where  the  tube  cross-section  goes  rapidly  to  zero  as  the  ray  approaches 
the  caustic,  and  ray  "C"  where  the  tube  cross-section  undergoes  rapid  expansion  as 
it  enters  the  shadow  zone  behind  the  caustic. 

RAY  TUBE  APPROACHING  CAUSTIC 


Results  for  wave  front  pressures,  p  , 

max 


the  caustic  are  shown  in  Figure  11. 
amplitude  factor, 


calculated  along  a  ray  tube  approaching 
These  results  are  presented  in  terms  of  the 


F  = 


r  A0(xh  I 

A  (x) 


(20) 


where  A(x)  »  ray  tube  cross-sectional  area  at  distance,  x,  along  the 

central  ray 

Aq(x)  »  ray  tube  cross-sectional  area  of  an  equivalent  non- 

refracted  ray  tube  (identical  Initial  boundary  rays)  at 
distance,  x,  from  the  charge,  i.e.,  in  this  case  the  area 

of  the  zone  of  a  sphere  of  radius,  x,  cut  by  the  boundary 
rays. 

(Text  continued  on  page  20)  H 
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FIG.  5  RAY  TUBE  AREA  AS  A  FUNCTION  OF  PATH  LENGTH 
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FIG,  6  RAY  TUBE  AREA  AS  A  FUNCTION  OF  PATH  LENGTH 
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PATH  LENGTH,  x  (FT) 

FIG.  7  RAY  TUBE  AREA  AS  A  FUNCTION  OF  PATH  LENGTH  —  7j  =  -15.7° 
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RAY  TUBE  AREA,  A(x)  (FT2) 
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RELATIVE  PATH  LENGTH 

FIG.  10  COMPARISON  OF  RAY  TUBE  AREA  FUNCTIONS 
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The  product 


F  *  'WlSO 


is  then  the  amplified  pressure  which  would  occur  along  the 
no  wave  overtaking,  i.e.,  in  the  acoustic  approximation— c 
ratio. 


ray  tube  if  there  were 

*  c  =  constant.  The 
o 


max 

F  x  ^Pmax^ISO 


therefore,  expresses  the  reduction  in  pmx  due  to  wave  overtaking  caused  by 
convergence  of  the  ray  tube  as  it  approaches  the  caustic. 


Since  we  are  interested  in  Isolating  the  effects  due  to  wave-overtaking,  both 
the  refracted  and  isovelocity  signatures  were  taken  from  calculations  which  are 
identical  except  for  the  cross-sectional  area  functions.  Thus,  in  the  comparison, 
small  errors  introduced  through  the  approximate  pressure-density  relation 
(equation  3)  and  the  shock  wave  simlltude  equations  (15)  and  (l6)  are  essentially 
cancelled . 

The  curves  shovn  in  Figure  11  were  calculated  for  the  ray,  7^  -15- 1°  (case 

B),  but  when  presented  in  this  manner  are  identical  to  those  for  the  ray, 

>.  *  -9.7°  (case  A).  The  three  curves  presented  in  Figure  11  were  calculated  for 

*  aT 

three  different  TMT  charge  weights — 56,  0.115,  and  2.6x10  lbs— which  yield 

spherical  wave  pressures  (no  refraction),  p  =175,  19,  and  0.46  psl,  at  256  ft 

max 

from  the  charge,  the  path  length  along  this  particular  ray  to  the  caustic.  The 

0.11$  and  2.6x10"^  lb  TOT  charge  weights  were  chosen  to  scale  the  caustic  location 
(x*2$6  ft  for  this  quarry  profile)  to  x=l80,000  ft  for  20  kiloton  and  900  lb  TNT 
charges,  respectively.  The  scaled  caustic  location  and  the  900-lb  charge  weight 
correspond  to  the  Sargasso  Sea  experiment  (Blatsteln,  1971). 


In  the  ocean  experiment,  the  caustic  was  a  convergence  zone  caustic.  When 
ws  did  these  calculations  we  had  hoped  that  the  results  would  be  approximately 
correct  for  rays  en  route  to  convergence  zone  caustics  in  the  ocean  as  well  as 
to  thermocllne -related  caustics.  But,  subsequent  calculations  using  ray  tubes  to 
convergence  zone  caustics  show  that  the  results  presented  in  Figure  11  are  not 
even  approximately  correct  for  convergence  zone  ray  tubes- -for  these,  the 
reduction  In  due  to  wave-overtaking  as  the  shock  wave  travels  to  the  caustic 

is  generally  much  greater,  sometimes  amounting  to  a  reduction  p _  of  25  or  30£. 

max 

We  also  point  out  that  the  results  presented  In  this  section  are  based  on  a  single 
experimentally  measured  sound  velocity  profile.  Thus,  the  results 
presented  In  Figure  11  may,  or  may  not,  be  typical  of  explosion  geometries  which 
form  thermocllne-related  caustics. 


Note  that  to  use  Figure  11  we  must  have  additional  knowledge  as  to  the 
refractive  effects  at  the  caustic.  In  practice  the  amplitude  factor,  F,  computed 
from  the  ray  tube  area  of  a  ray  theory  computation  is  of  no  use  in  estimating  the 
effects  of  wave-overtaking  using  Figure  11,  since  In  the  area  of  interest  wave 
diffraction  has  already  negated  such  a  computation.  Consequently,  to  use  Figure  11 
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one  must  resort  to  experimental  data  or  a  modified  ray  theory  computation  to 
obtain  a  value  for  F,  using  the  equation 


F  2: 


p  max^max^ISO 


(21) 


where  p'  is  some  value  of  p  from  which  one  can  estimate  an  effective  ray 
e  max  max 

tube  area  and  (p  is  the  pressure  of  a  spherical  wave  which  has  traveled  the 

maxISO 

same  path  length. 


In  the  quarry  and  ocean  experiments  to  date  increases  in  peak  pressure  of, 

say,  4,  6,  or  10— but  no  greater  than  a  factor  of  10 — have  been  observed  at 

caustics.  This  Is  also  true  of  the  results  of  modified  ray  theory  computations 

based  on  sound  velocity  profiles  taken  from  these  experiments.  Thus,  according  to 

the  curves  shown  in  Figure  11  reductions  in  p  due  to  wave-overtaking  of  up  to 

max 

4  or  5 i>  occur  as  the  shock  wave  approaches  this  particular  caustic. 

RAY  TUBE  EXPANDING  IN  SHADOW  ZONE 


Figure  12  shows  the  calculated  pressure  as  a  function  of  distance  along  the 

ray,  7^=  -15.7°  (case  C),  and  compares  it  with  the  pressure  calculated  for  a 

spherically  expanding  6hock  wave.  For  this  case  the  deviation  in  p  from 

max 

spherical  wave  pressure  can  be  predicted  quite  well  by 


max 


x  max^ISO 


(22) 


which  is  also  shown  in  Figure  12  (dashed  line).  Out  to  350  ft  equation  22  gives 

p  to  approximately  ±l£--at  x=600  ft  equation  22  predicts  a  value  of  p  about 
max  max 

4$  lover  than  computed  by  the  finite  amplitude  calculation. 


5.  DISCUSSION  AND  CONCLUSIONS 


As  pointed  out  in  the  introduction,  existing  solutions  to  the  refraction 
problem  (ray  theory  and  modified  ray  theory)  do  not  take  into  account  wave- 
overtaklng  and  molecular  absorption.  However,  the  computed  pressures  are 
generally  corrected  afterwards  in  the  same  ratio  as  the  correction  for  wave¬ 
overtaking  and  absorption  estimated  for  a  spherically  diverging  wave  which  has 
traveled  the  same  distance.  The  purpose  of  this  paper  vms  to  check  this 
approximate  method  of  accounting  for  wave-overtaking.  To  do  this  we  made  an 
essentially  exact  calculation  Including  the  effects  of  wave-overtaking  (but  not 
absorption)  for  an  approximate  problem:  namely,  the  propagation  of  an  underwater 
explosion  shock  wave  along  a  duct  formed  by  two  adjacent  rays  calculated  by  ray 
theory. 


The  computational  results  (Section  4)  are  for  two  extreme  cases— a  converging 
ray  tube  approaching  a  caustic  and  one  undergoing  rapid  expansion  as  it  enters 
the  shadow  region  behind  the  caustic.  These  were  selected  from  a  particular 
sound  velocity  profile  forming  a  thermocllne -related  caustic.  For  these  cases 
the  approximate  procedure  (estimating  the  wave-overtaking  effects  from  the 
spherical  wave)  Introduced  a  maximum  error  of  ±5$,  But,  as  pointed  out  in 
Section  4,  this  result  may  or  may  not  be  typical  of  ray  tubes  to  thermooline-related 
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PATH  LENGTH,  x  (FT) 

FIG.  12  CALCULATED  .  JITIAL  PRESSURE,  p  ,  ALONG  AN  EXPANDING  RAY  TUBE 

max 


caustics.  Also,  subsequent  calculations  covering  a  range  of  sound  velocity 

profiles  yielding  convergence  zone  caustics  resulted  In  considerably  greater 

reductions  In  p  due  to  wave-overtaking.  For  these  cases  most  of  the 
max 

attenuation  ocurred  along  the  portion  of  the  ray  path  in  the  deep  ocean  velocity 
gradient. 

Figure  11  shovs  that  for  ray  tubes  approaching  this  particular  thermocline- 
related  caustic  there  is  increased  attenuation  of  due  to  vave-overtaklng  as 

the  pressure  builds  up.  This  additional  attenuation  is  dependent  on  the  size  of 
the  charge- -or  the  pressure  level  of  the  wave  approaching  the  caustic--in 
addition  to  the  variations  in  ray-tube  cross-sectional  area  between  the  charge 
and  the  caustic.  The  moderate  magnitude  of  the  vave-overtaklng  corrections  to 

p  would  appear  to  Justify  the  current  practice  of  predicting  explosion 

max 

pressures  at  and  near  caustics  with  modified  ray  theory  computations.  Guch 
results  can  then  be  adjusted,  if  necessary,  for  attenuation  due  to  wave-overtaking 
occurring  along  the  path  to  the  caustic  by  the  method  presented  in  this  report. 

Several  results  from  subsequent  calculations  of  wave  propagation  along  ray 
tubes  to  convergence  zone  caustics  in  the  ocean  have  been  mentioned  above. 

Further  results  and  details  of  these  calculations  will  be  given  in  a  subsequent 
report. 
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